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Abstract 

Brucellosis is a worldwide distributed zoonosis that causes important economic losses to animal production. In Brazil, 
information on the distribution of biovars and genotypes of Brucella spp. is scarce or unavailable. This study aimed (/') to 
biotype and genotype 137 Brazilian cattle isolates (from 1977 to 2008) of B. abortus and (;'/) to analyze their distribution. 6. 
abortus biovars 1, 2 and 3 (subgroup 3b) were confirmed and biovars 4 and 6 were first described in Brazil. Genotyping by 
the panel 1 revealed two groups, one clustering around genotype 40 and another around genotype 28. Panels 2A and 2B 
disclosed a high diversity among Brazilian B. abortus strains. Eighty-nine genotypes were found by MLVA16. MLVA16 panel 1 
and 2 showed geographic clustering of some genotypes. Biotyping and MLVA16 genotyping of Brazilian B. abortus isolates 
were useful to better understand the epidemiology of bovine brucellosis in the region. 
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Introduction 

Brucellosis is a worldwide distributed zoonosis, responsible for 
considerable economic losses due to abortion and culling of 
infected animals [1]. The genus Brucella is currently composed of 
ten species preferentially associated with different host: B. abortus 
(cattie), B. canis (canids), B. ceti (cetaceans), B. melitensis (sheep and 
goats), B. microti (Microtus arvalis), B. neotomae (Motoma lepida), B. ovis 
(sheep), B. pinipedialis (pinnipeds), B. suis (pigs), and B. inopinata 
(man?) [2]. 

Bacteriological culture is considered the gold standard for 
brucellosis diagnosis. Due to its high specificity and ability to 
identify different species and biovars (bv), bacteriological culture is 
very important for epidemiological studies, allowing the determi- 
nation of the way the disease spreads in a region and the possible 
sources of infection [3,4]. 

Brucella abortus biovars have distinct geographic distributions [5] 
Brucella abortus bvl and B. abortus bv2 are worldwide distributed, 
while B. abortus bv3 is found predominantly in Italy, India, Egypt, 
and Africa. B. abortus bv5 is most commonly found in Germany 
and United Kingdom [5], but has also been observed in France 
[6] . B. abortus biovars 4 and 6 have been reported on Mexico and 
France, though less frequendy than bvl, 2, and 3 [6,7], B. abortus 



bv4 was also reported in Canada [8], Chile, Ecuador, and Cuba 
[9]. In India, B. abortus bvl is the most frequent, followed by B. 
abortus bv3, although bv4, 6, and 9 were also found [10]. 

There have been just few studies concerning the identification of 
the biovars of Brucella spp. in Brazil describing B. abortus bvl, 2 and 
3 and B. suis bvl [9,11,12]. Also, B. canis and B. ovis were found 
infecting domestic animals [3] . Brucella melitensis was not reported 
in Brazil [3,4]. 

Although the great homogeneity of the genus Brucella [13], 
conserved loci containing repeated sequences were found. Analysis 
of the number of repeats of these loci (Variable Number of 
Tandem Repeats - VNTR) is used as tool in the study of 
brucellosis [14-17]. Also, VNTR analysis of various loci (Multiple 
Locus VNTR Analysis - MLVA) is very important for investiga- 
tions of infection by Brucella spp. and was applied to epidemio- 
logical control and surveillance of brucellosis [14,17-22]. More- 
over, this assay proved to be a valuable tool to determine 
relationships among Brucella isolates from different animal species 
and from humans, as well as for epidemiological trace-back studies 
[23,24]. 

Thus, the aims of the present study were (i) to biotype and 
genotype Brazilian cattle B. abortus strains and (a) to analyze their 
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distribution to support the Programa National de Controle e 
Erradicacao de Brucelose e Tuberculose (PNCEBT) [25]. 

Materials and Methods 

Bacterial Strains 

One hundred and thirty-seven Brazilian cattle B. abortus isolates, 
obtained from 1977 to 2008, were used. They were from the 
collections of Laboratorio de Bacteriologia Aplicada, Escola de 
Veterinaria, UFMG; Laboratorio National Agropecuario, Insti- 
tute de Pesquisas Veterinarias Desiderio Finamor, and Instituto 
Biologico de Sao Paulo. These isolates were from Minas Gerais 
(38), Rio Grande do Sul (31), Para (31), Sao Paulo (15), Tocantins 
(14) and Santa Catarina (8) (Figure 1). These strains represent the 
great majority of B. abortus isolates currently available in Brazil. 

Reference strains, that were used as control in different 
procedures, were: B. abortus biovar (bv) 1 544 = ATCC 23448 T , 
B. abortus bv2 ATCC 23449, B. abortus bv3 Tulya = ATCC 23450, 
B. abortus bv4 292 = ATCC 23451, B. abortus bv5 3 196 = ATCC 
23452, B. abortus bv6 870 = ATCC 23453 and B. abortus bv9 
C68=ATCC 23455, B. abortus bvl B19 USDA, B. abortus bvl 
RB51, B. melitensis bvl 16M = ATCC 23456 T , B. ovh Reo 198, B. 
suis bvl 1 330 = ATCC 23444, Escherichia coli ATCC 25922, E. coli 
B41, Listeria monocytogenes ATCC 19115, Pseudomonas aeroginosa 
ATCC 27853, Salmonella enterica biovar Typhimurium ATCC 
14028 and Staphylococcus aureus ATCC 29213. 



Identification and Biotyping of B. abortus 

All strains were identified by routine test [26,27]. Isolates were 
also tested by bcsp31 -PCR [28], AMOS-PCR [28], AMOS-ERY- 
PCR [30] for classification of subgroups 3a and 3b, and omp2b- 
RFLP-PCR with restriction by Taql to confirm the identity of B. 
abortus bv4 strains [31]. 

MLVA16 Genotyping 

DNA from each strain was submitted to MLVA typing, using a 
subset of 16 tandem repeat loci (MLVA16) [15,17]. 

Genotype Analysis 

Band size estimates were converted into number of repeat units 
for each locus [15], with the software BioNumerics 6.1 (Applied 
Maths, Sint-Martens-Latem, Belgium). Clustering analysis also 
was performed with the software BioNumerics 6.1 based on 
UPGMA [15]. The diversity index HGDI was calculated [32]. 
Genotypes obtained were compared to those deposited in the 
MLVAbank for Bacterial Genotyping (http://mlva.u-psud.fr/ 
mlvav4/genotyping/index.php). The minimum spanning tree 
presented is the one with the highest overall reliability score and 
was calculated using UPGMA associated with the priority rule and 
the bootstrap resampling (BioNumerics 6.1). 
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Figure 1. Sampling and distribution of biovars and MLVA16 Panel 2 genotypes of B. abortus isolates from Brazil, 1977 to 2008. 

doi:10.1371/journal.pone.0081152.g001 
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Results 

Biotyping and Identification of B. abortus Isolates 

All the 137 isolates were confirmed to be Brucella spp. by bcsp'31- 
PCR and to be B. abortus by biochemical tests, AMOS and 
AMOS-ERY-PCRs. Ninety-two (67.1%) of the 137 isolates were 
identified as B. abortus bvl, B. abortus bv2 or B. abortus bv4 by 
AMOS-PCR. These were found in all studied States, being bvl 
most frequent in Rio Grande do Sul (26/72), followed by Minas 
Gerais (23/72), Sao Paulo (12/72), Santa Catarina (6/72), Para 
(3/72) and Tocantins (2/72). B. abortus bv2 was most frequent in 
Minas Gerais (10/16), followed by Santa Catarina (2/16), Rio 
Grande do Sul (2/16), Para (1/16) and Tocantins (1/16), while 
strains of bv3 were mainly observed in Para (16/28), followed by 
Tocantins (9/28), Minas Gerais (2/28) and Sao Paulo (1/28). 
Similarly, bv6 strains were also more common in Para (11/20), 
followed by Minas Gerais (3/20), Rio Grande do Sul (2/20), Sao 
Paulo (2/20) and Tocantins (2/20). Interestingly, only one strain of 
B. abortus bv4 was observed in Santa Catarina state. AMOS-PCR- 
negative strains were found in all States except Santa Catarina and 
most frequently in Para and Tocantins. 

B. abortus biovars found by biochemical testing were 1,2, 3, 4, 
and 6. The most frequent was B. abortus bvl, followed by bv3, bv6, 
and bv2; the least frequent was bv4. The HGDI calculated for B. 
abortus biovar identification was 0.65. All B. abortus bv3 strains were 
classified as subgroup 3b by AMOS-ERY-PCR. 

The unique B. abortus bv4 isolate, strain 16/02, isolated in Rio 
Grande Sul, was confirmed by om/>2b-RFLP-PCR [31]. 

MLVA16 Typing and Clustering of B. abortus Isolates 

MLVA16 panel 1 identified nine genotypes (Table 1), four (28, 
32, 33 and 40) previously described [15] and five new ones (I to V). 
These genotypes had one or more loci of panel 1 minisatellites with 
a number of repeat units different from reported ones [15]. Among 
genotypes observed for panel 1 genotyping the most frequent was 
genotype 28 found in all studied states, being more common in 
Minas Gerais (33), followed by Rio Grande do Sul (27/82), Sao 
Paulo (9/82), Santa Catarina (8/82), Tocantins (3/82) and Para 
(2/82). The genotype 40 was more observed in Para (26/45), 
followed by Tocantins (9/45), Minas Gerais (5/45), Sao Paulo (3/ 
45) and Rio Grande do Sul (2/45). Whereas, the three strains of 
genotype 33 were observed in Tocantins (2/3) and Sao Paulo (1/3) 
and genotype 32 was found only in Rio Grande do Sul (1/1). 
Regarding new genotypes of panel 1, the profile I and II were 
observed only in Tocantins state, while patterns III and IV were 
found respectively in Para and Rio Grande do Sul, all with one 
representative each. With two representative, genotype V was 
observed only in Sao Paulo state. The HGDI for panel 1 
genotyping was 0.54. 

Panel 2A identified 13 genotypes, five (8, 30, 32, 33 and 34) 
previously described [15] and eight new ones (2A1 to 2A8). Panel 
2A HGDI was 0.74. Genotype 33 was the most frequent (62/137), 
followed by 32 (28/137), 34 (13/137), 8 (10/137) and 30 (8/137). 
Among newly described genotypes, 2A3 was the most frequent (4/ 
137), followed by 2A2 and 2A8 (3/137). 

In panel 2B, 63 different genotypes were found, 13 of which had 
been previously described (18, 29, 30, 45, 52, 53, 58, 64, 66, 87, 
88, 90 and 91) and 50 new genotypes, identified as 2B1 to 2B50. 
The HGDI for these loci was 0.97. Among these genotypes, the 
most frequent was genotype 91 (1 1/137), followed by 19 (9/137), 4 
(8/137), 5 (6/137), 18 and 8 (5/137). The genotypes 30, 53, 17, 
22, and 48 were found each in four of the 137 strains and 
genotypes 52, 58, 90, 16, 30, 44 and 46 in three of them. 



Table 2 shows HGDI for the 16 loci based on biovars of studied 
B. abortus strains. Panel 1 had HGDIs below 0.16. Panel 2A loci 
with higher diversity were Bruce 19 and Bruce21. HGDIs for B. 
abortus bv2 (0.74) and bv6 (0.73) were the highest in panel 2A. 
Panel 2B had the highest HGDI among all biovars (0.93 to 0.96); 
with B. abortus bvl presenting the highest one. HGDI for MLVA16 
(panels 1, 2A and 2B) was 0.98, differentiating 89 genotypes 
among the 137 studied strains (A to CY). 

The generated Minimum Spanning Tree (MST), using cate- 
gorical and parsimonious coefficients, clearly separated the two 
groups (Figure 2). 

Among MLVA16 genotypes found, the most frequent was 
genotype J (9/137; 6.6%), found in Minas Gerais, Sao Paulo, and 
Santa Catarina. Other frequent genotypes were D (5.8%; 8/137), 
found in Minas Gerais and Rio Grande do Sul; AB (3.6%; 5/137) 
and U (2.9%; 4/ 1 37), both found in Para and Tocantins; H and R 
(2.2%; 3/137) from Minas Gerais, Para and Tocantins, and Minas 
Gerais and Sao Paulo, respectively. Genotypes T, AC and BE 
(1.5%; 2/137) were found in Para and Tocantins, and genotypes L 
and M (2.2%) in Minas Gerais and Santa Catarina, and in Rio 
Grande do Sul and Sao Paulo, respectively. Table SI summarize 
the genotypic profile in MLVA16, genotypes in panel 1, alleles in 
each MLVA16 loci, biovar and field information of 137 cattle B. 
abortus isolates. 

All new genotypes found in this study were included in 
MLVAbank (http:/ / mlva.u-psud.fr/ mlvav4/ genotyping/index.php) 
in a cooperative database called Brucella_Brazil. 

Comparison among Typing Methods 

Comparison of biochemical test results with panel 1 genotyping 
and AMOS-PCR showed compatible results by the three 
techniques in 68.6% (94/137) of the cases (Table 1). Twenty- 
three strains (16.8%) biochemically classified as B. abortus bv3 were 
from genotype 40 and negative by AMOS-PCR. Three strains 
biochemically classified as B. abortus bv3 were positive by AMOS- 
PCR and also from genotype 40. One strain phenotypically 
identified as B. abortus bv4 was AMOS-PCR positive and from 
genotype 32. Five strains (10.9%) that were AMOS-PCR positive 
B. abortus bv2 were typed as genotype 28. Another AMOS-PCR 
positive B. abortus bv2 was from genotype 33. 

Discussion 

Brucella abortus biovars found in this study were 1, 2, 3, 4, and 6 
(Table 1), the latter two being reported for the first time in Brazil. 
B. abortus bvl was the most frequent biovar isolated, as also 
described by previous studies in Brazil [3,12,33] and most regions 
of the world [4, 1 8] . Interestingly, our B. abortus bv2 isolates did not 
need addition of serum to culture medium for their isolation, as it 
is usually described for strains of this biovar [27]. This 
phenomenon was also observed for other B. abortus bv2 strains 
from Brazil [12]. Of the two biovars that were described for the 
first time in Brazil (bv4 and bv6), only one strain of B. abortus bv4 
was isolated from Rio Grande do Sul. Conversely, B. abortus bv6 
was observed in all the States sampled except Santa Catarina. 

In our study, the high isolation frequency of A abortus bv 1 and 
bv2 could be explained by the continuous importation of cattle, 
since the beginning of colonization [34] , from regions where these 
two biovars are or were endemic [ 1 ,4, 1 0] . As all isolates of B. 
abortus bv3 were from sub-group 3b, they were more closely related 
to B. abortus bv3 of European than of African origin [30,35], which 
is in accordance to the origin of most cattie imported into Brazil, 
that came from Europe and India [10,34].The unique B. abortus 
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Table 1. Biotyping, AMOS-PCR and genotyping based on panel 1 of MLVA16 of the strains of B. abortus isolated from cattle in 
Brazil, 1977 to 2008. 





Biovar 


AMOS-PCR 


Genotype 3 


















28 (%) 


32 (%) 33 (%) 


40 (%) 


1 (%) 


II (%) 


III (%) 


IV (%) 


V (%) 


Total (%) 


1 


Pos 


67 (48.9) 


2 (1.5) 


- 


- 


- 


- 


1 (0.7) 


2(1.5) 


72 (52.6) 




Neg 


b 


- 


- 


- 


- 


- 


- 


- 


0 (0.0) 


2 


Pos 


15 (10.9) 


1 (0.7) 


- 


- 


- 


- 


- 


- 


16 (11.7) 




Neg 


















0 (0.0) 


3 


Pos 






3 (2.2) c 












3 (2.2) 




Neg 






23 (16.8) 




1 (0.7) 


1 (0.7) 






25 (18.6) 


4 


Pos 




1 (0.7) 














1 (0.7) 




Neg 


















0 (0.0) 


6 


Pos 


















0 (0.0) 




Neg 






19 (13.9) 


1 (0.7) 










20 (14.6) 


Total 




82 (59.9) 


1 (0.7) 3 (2.2) 


45 (32.8) 


1 (0.7) 


1 (0.7) 


1 (0.7) 


1 (0.7) 


2 (1.5) 


137 (100.0) 



a Panel 1 of MLVA16 is based on loci Bruce06, Bruce08, Brucell, Bruce12, Bruce42, Bruce43 and Bruce55 [15]. 
b n = 0. 

c result not expected. 

doi:1 0.1 371 /journal.pone.0081 1 52.t001 



bv4 isolate were from Rio Grande do Sul, at the border with B. abortus bv6 was not previously reported in South America. 

Argentina, where this biovar has already been described [4]. The highest frequency of this biovar was found in Para, where 

55% of the isolates were of B. abortus bv6. Although the way this 



Table 2. Diversity index (HGDI) calculated for each locus and for panels 1, 2A, 2B and MLVA16 according to the biovars of 8. 
abortus isolated from cattle, Brazil, 1977 to 2008. 



Group/Locus 


Biovar 1 




Biovar 2 




Biovar 3 




Biovar 6 




Types" 


HGDl" 


Types 


HGDI 


Types 


HGDI 


Types 


HGDI 


Panel 1 


5 


0.16 


2 


0.07 


3 


0.14 


2 


0.1 


Bruce06 


2 


0.05 


2 


0.12 


2 


0.07 


2 


0.09 


Bruce08 


2 


0.03 


1 


0.00 


2 


0.07 


1 


0.00 


Brucel 1 


2 


0.03 


1 


0.00 


1 


0.00 


2 


0.09 


Brucel2 


3 


0.08 


1 


0.00 


1 


0.00 


2 


0.09 


Bruce42 


1 


0.00 


1 


0.00 


1 


0.00 


1 


0.00 


Bruce43 


1 


0.00 


1 


0.00 


1 


0.00 


1 


0.00 


Bruce45 


2 


0.03 


1 


0.00 


1 


0.00 


2 


0.09 


Bruce55 


1 


0.00 


1 


0.00 


2 


0.07 


1 


0.00 


Panel 2A 


7 


0.42 


4 


0.74 


4 


0.62 


7 


0.73 


Brucel 8 


3 


0.18 


1 


0.00 


1 


0.00 


4 


0.27 


Brucel 9 


2 


0.13 


2 


0.50 


3 


0.37 


3 


0.40 


Bruce21 


2 


0.26 


3 


0.62 


2 


0.35 


5 


0.54 


Panel 2B 


38 


0.96 


10 


0.93 


18 


0.95 


12 


0.95 


Bruce04 


4 


0.49 


2 


0.50 


9 


0.88 


8 


0.87 


Bruce07 


3 


0.46 


2 


0.12 


2 


0.07 


3 


0.33 


Bruce09 


4 


0.08 


1 


0.00 


2 


0.07 


2 


0.09 


Brucel 6 


6 


0.70 


3 


0.59 


2 


0.35 


2 


0.18 


Bruce30 


6 


0.62 


4 


0.67 


3 


0.20 


2 


0.18 


MLVA16 


48 


0.93 


14 


0.98 


21 


0.97 


16 


0.98 



a Number of genotypes; 

b Diversity index calculated according to Hunter and Gaston (1988) [32]. 
doi:1 0.1 371 /journal.pone.0081 1 52.t002 
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Minas Gerals - MG 

1 1 

| Sao Paulo ■ SP 

| Para - PA 

Tocantins - TO 

Rio Grande do sul - RS 

{^P} Santa catarlna ■ SC 

8. abonus reterence strains 

Figure 2. Minimum Spanning Tree (MST) analysis of B. abortus isolates using the MLVA1 6 data. The minimum spanning tree presented is 
the one with the highest overall reliability score and was calculated using UPGMA associated with the priority rule and the bootstrap resampling for 
strains of Brucella abortus isolated from cattle in Brazil, 1977 to 2008, the reference strains for each biovar of S. abortus and the vaccine strains S19 and 
RB51. Numbers in each clonal complex represent the genotype on Panel 1. Information on the origin of the isolates was color labeled. 
doi:10.1371/journal.pone.0081152.g002 



biovar was introduced into Brazil is not clear, B. abortus bv6 is one 
of the most frequent biovars of B. abortus in cattle in India, which is 
origin of the Zebu breeds that predominate in Para and most of 
the country [10,34]. 

Although MLVA16 could correcdy group the strains by species 
of Brucella, it was not able to identify Brucella spp. biovars, probably 
due to the high rates of mutation at some loci. Therefore, this 
technique is not an adequate substitute for classic biovar 
identification [15], which is supported by this study, as three 
genotypes (28, 33 and 40) were found in B. abortus from more than 
one biovar (Table 1). 

The most frequent genotype among the studied Brazilian B. 
abortus isolates was genotype 28, found in all six States. In Santa 
Catarina, it was the only genotype found among the eight isolates 
of B. abortus. Genotype 28, generally reported to be associated with 
strains of B. abortus bvl and 4 from USA, Costa Rica, Switzerland, 
France, Italy, and Germany, was observed in strains of B. abortus 
bvl and in almost all (15/16) of B. abortus bv2 strains examined. 

Genotype 40, first reported from strains of B. abortus bv6 from 
Africa, was also very frequently found in Brazil and is the most 
common among strains of A abortus bv6 (95.0%) and B. abortus bv3 
(92.9%). This explains the geographic distribution of this 



genotype, as these two biovars were also the most common in 
Para and Tocantins. 

The unique B. abortus bv4 strain, isolated from Rio Grande do 
Sul, was from genotype 32, previously reported among German B. 
abortus bvl strains. Genotype 33 strains, described in B. abortus bvl 
strains from France, were observed in this study in B. abortus bv3 
strains from Para (2) and Sao Paulo (1). The lack of sound 
epidemiological data on those isolates precludes any deeper 
inference on the way this genotype circulates in Brazil. 

Among the newly described panel 1 genotypes, the major 
differences of genotype I are the presence of 1 1 repeat units at locus 
Bruce 12 and six repeat units at locus Bruce45, whereas genotype II 
had only one repeat unit at locus Bruce55. Genotype III had as 
leading feature the presence of just four repeat units at locus 
Bruce08 and genotype IV, observed in a B. abortus bvl isolated 
from Rio Grande do Sul, was mainly characterized by four repeat 
units at locus Bruce45. Genotype V strains are characterized by the 
presence of 13 repeat units at locus Bruce 12. All of these new 
genotypes were still represented by just one strain each, exception 
for genotype V, which was observed in two strains from Sao Paulo. 
Additional systematic studies in the regions where those genotypes 
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were observed could help to understand their epidemiological 
importance. 

The HGDI calculated for MLVA16 panel 1 (0.54) was lower 
than that calculated for biotyping (0.65), indicating a lower 
discriminatory efficacy of MLVA16 panel 1, due to the presence of 
more than one biovar in the same genotype by this panel (Table 1). 
The HGDI calculated for MLVA16 panel 2A (0.74) is similar to 
that found by Al Dahouk et al. (2007) [17], who analyzed 128 
human isolates of B. melitensis and found that the mean diversity 
index in panel 2A was less than 0.75. This panel has mainly helped 
in the differentiation of strains of bv2 and 6 (Table 2). The high 
HGDI calculated for MLVA16 panel 2B, which has the most 
variable loci, with diversity indexes generally higher than 0.80, 
demonstrates that it can be used, together with MLVA16 panel 
2A, for identifying strains associated with disease outbreaks 
[14,15,17]. Comparing HGDI for MLVA16 panels 2A and 2B 
calculated in the present study and those described in the literature 
[15,17,19-21] leads to the conclusion that these loci were effective 
for differentiating the Brazilian B. abortus isolates, especially when 
they are compared to biotyping. 

Genotyping based on MLVA16 identified 89 previously unde- 
scribed genotypes (http://mlva.u-psud.fr/mlvav4/genotyping/ 
index.php) among Brazilian B. abortus strains. The difference 
between the number of genotypes identified by MLVA16 panel 1 
and by MLVA16 panel 2 is a consequence of the hypervariability of 
the loci in this second panel [15,17]. 

The results of genotyping by MLVA16 show that the same 
genotypes tend to circulate in neighboring States, such as 
genotypes AB, U, T, AC, and BE in Para and Tocantins in 
North Region, and genotypes J, D, L, and M in Minas Gerais, Sao 
Paulo, Santa Catarina, and Rio Grande do Sul in Southeast and 
South regions. However, some genotypes, as genotype H, found in 
North (Para and Tocantins) and Southeast (Minas Gerais) regions, 
were present in distant and unrelated regions, thus requiring 
newer studies to clarify their epidemiology. 

The MST (Figure 2) showed that the B. abortus isolates from six 
Brazilian States were displayed into two well-defined clusters, with 
51.7% similarity. Cluster I included B. abortus strains of genotypes 
40, 33 (2/3), II and III, distributed in five of the six States, except 
Santa Catarina. Cluster I also included reference strains of B. 
abortus bv5 3196, B. abortus bv6 870 and B. abortus bv9 C68. Cluster 
II included B. abortus strains from genotypes 28, IV, 33 (1/3) and 
V, as well as reference strains of B. abortus bvl 544, B. abortus bv2 
86/8/59 and B. abortus bv4 292 and vaccine strains B. abortus bvl 
S19 and B. abortus bvl RB51. Clustering was performed basically 
from locus Bruce04, which was the most variable, followed by locus 
Bruce 19 and then by locus Bruce30, first grouping strains with 
higher number of repeat units and then those that had fewer 
repeat units. Similarly, MLVA analysis of the Spanish human B. 
melitensis strains showed a high diversity of genotypes with locus 
Bruce04 featuring the highest HGDI [20] . 

Two strains appeared apart from the main clusters: B. abortus 
bv3 reference strain Tulya, isolated from Africa, which is 
MLVA16 genotype 207, and B. abortus bv6 strain 95, isolated 
from Tocantins, which is MLVA 16 genotype AR. This shows that 
these strains are more distantly related to the other Brazilian B. 
abortus isolates (Figure SI). Likewise, this can be observed in MST 
(Figure 2), which also shows that strain Tulya is distantly linked to 
cluster II (gathered around genotype 28) and that strain 95 is 
distandy linked to cluster I (gathered around genotype 40). 

Typing the strains by biotyping, genotyping by MLVA 16 panel 
1 and AMOS-PCR gave compatible result in the majority of the 



cases. Most of the genotypes that were observed in the present 
study and were not previously described in a particular biovar 
(Table 1) are probably due to small numbers of isolates from these 
biovars that were already analyzed by MLVA 16, being the present 
ones new descriptions of these genotypes in different biovars. 

B. abortus bv3 is a heterogeneous biovar, divided into subgroups 
3a and 3b [15,30] and the MST (Figure 2) depicts these 
differences. All Brazilian B. abortus bv3 strains were closer to B. 
abortus bv6 reference strain 870 than to the B. abortus bv3 reference 
strain Tulya, which did not cluster with other strains. These results 
agreed to published reports, which showed that strains of B. abortus 
bv3 subgroup 3b present a pattern similar to that of strains of A 
abortus bv5, B. abortus bv6 and B. abortus bv9 [30]. 

Interestingly, three strains of B. abortus bv3, genotype 40, 
isolated from Para, showed positive results in AMOS-PCR, which 
was described as able of detecting only biovars 1, 2, and 4 of B. 
abortus [29]. However, occasionally, atypical isolates of Brucella 
spp., that do not match identification tables and keys [36], has 
been described. Banai et al. (1990) [37] found B. mellitensis bvl 
strains that grew with acid fuchsin and thionine in two regions of 
Israel, suggesting evolution of a new variant. Corbel (1991) [38] 
also found isolates of A mellintensis in Italy, Kuwait, Saudi Arabia, 
Zimbabwe, India, and Germany that were sensitive to thionine 
(20 ng/mL). Garcia et al. (1988) [39] isolated atypical strains of B. 
abortus from seven infected cattie herds in Ontario (Canada) that 
were sensitive to thionine and fuchsin, similar to B. abortus bv2, but 
they agglutinated in anti-M serum, which is characteristic of B. 
abortus bv4. 

When these atypical isolates appear, it is necessary to determine 
if their appearance is a reproducible phenomenon, indicating a 
local taxonomic pattern, and not a laboratory error. If reproduc- 
ibility is proven, this can be used as an epidemiological marker 
[27]. Bricker (2002) [36] reported that due to the small number of 
differences among species and biovars of Brucella spp., the smallest 
mutation could result in conflicting data, making interpretation of 
the characterization of the isolates difficult. The three strains of B. 
abortus bv3, genotype 40, that were reactive in AMOS-PCR were 
isolated from the municipality of Conceicao do Araguaia, Para, 
within an interval of four years between the isolation of the first 
strain and the other two; which were isolated within an interval of 
five months. This could be an indication that these markers are 
characteristic of atypical strains from this region. 

Biotyping and MLVA 1 6 genotyping of Brazilian isolates of B. 
abortus will allow the construction of a database of genotypes of 
Brucella spp. from Brazil, and, in the future, for Latin America, 
which certainly will contribute to a better understanding of the 
epidemiology and control of bovine brucellosis in the region. 
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